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Abstract — High-statistics data on near-horizontal muons 
collected with Russian-Italian coordinate detector DECOR 
are analyzed. Precise measurements of muon angular distri- 
butions in zenith angle interval from 60° to 90° have been 
performed. In total, more than 20 million muons are se- 
lected. Dependences of the absolute integral muon intensity 
on zenith angle for several threshold energies ranging from 
1.7 GeV to 7.2 GeV are derived. Results for this region 
of zenith angles and threshold energies have been obtained 
for the first time. The dependence of integral intensity on 
zenith angle and threshold energy is well fitted by a simple 
analytical formula. 



I. Introduction 

STUDIES of angular and energy dependence of muon 
flux at the Earth's surface give important information 
as about processes of muon generation and propagation in 
the atmosphere so about primary cosmic rays. Measure- 
ments of muon flux at large zenith angles up to 90° are 
especially actual since primary particles for such muons 
have higher energies than in the vertical direction. Ex- 
perimental studies of muon intensity at large zenith an- 
gles at the ground level can be conditionally separated in 
two groups: measurements of muon integral intensity with 
threshold energies less than 1 GeV [1]— [8] and investiga- 
tions of integral and differential muon spectra for muon 
energies higher than 10 GeV (see review [9]). Regions of 
measurements of muon spectrum at large zenith angles are 
presented in Fig. 1. It is remarkable that for threshold en- 
ergies from 1 GeV to 10 GeV and zenith angles 60°<6><90° 
muon intensity data are absent. 

To explore this region, a setup capable to measure near- 
horizontal muon flux at different threshold energies with a 
good angular accuracy of track reconstruction is needed. 
Coordinate detector DECOR, which is a part of experi- 
mental complex NEVOD situated in MEPhI (Moscow), is 
such a detector. Regions of threshold energies and zenith 
angles accessible for DECOR and analyzed in this work 
are shown by the dashed areas in Fig. 1. 
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Fig. 1 

Regions of muon spectrum measurements at large zenith 
angles. Symbols represent the measurements of integral 
intensity; shaded areas are the regions of differential 
spectrum measurements. dashed areas are the regions 
investigated in this work. 



II. Experimental setup and events selection 

Experimental complex NEVOD includes a water 
Cherenkov calorimeter NEVOD [10] with sensitive volume 
2000 m 3 equipped with quasispherical modules of PMTs, 
and large-area (~ 110 m 2 ) coordinate detector DECOR 
[11] (Fig. 2). Eight supermodules (SM) of DECOR are 
situated in the gallery around the water tank, and four 
SM on its cover. SM of side part of DECOR represents 
eight parallel planes with sensitive area 3.1 m x 2.7 m, 
suspended vertically with 6 cm distance from each other. 
These planes consist of 16 chambers which contain 16 tubes 
with inner cross-section 0.9 cm x 0.9 cm. Chambers are 
operated in a limited streamer mode and are equipped with 
two-coordinate external strip read-out system. Thus, co- 
ordinates of passing particle can be obtained for each plane 
with spatial accuracy of muon track location ~ 1 cm. First 
level trigger is formed when there are at least two even and 
two odd triggered planes in a given SM. 

For the analysis, particles passing through two SM situ- 
ated at different sides of the water pool were selected. Dif- 
ferent pairs of SM correspond to different values of thresh- 
old energy. Accuracy of zenith angle reconstruction for 
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Fig. 2 

Experimental complex NEVOD-DECOR. 



tracks passing through selected SM pairs is 0.3° —0.5°. 
Selection procedure includes the following conditions. 

• "OneTrack" criterion: two tracks reconstructed from 
data of different supermodules must coincide within 
5° cone. In this case the tracks in separate SMs are 
considered as tracks of the same particle. Straight 
line connecting the middles of two reconstructed track 
segments is taken as the trajectory of the particle. 

• The events in which muon passed closer than 3 cm 
from the boundary of SM are rejected in order to de- 
crease the edge effects. 

« There must be two and only two track projections 
(X,Y) in each SM for unambiguous reconstruction of 
geometrical characteristics of muon track (the absence 
of accompanying particles). 
Data collected over a period from December 2002 to 
June 2003 are analyzed. Total time of registration is equal 
to 3390 hours. The total number of selected events is more 
than 20 millions. 

III. Results 

Threshold energy E m i n of muons passing through se- 
lected pair of SM is calculated by means of range-energy 
tables [12] . It is calculated for each selected event, and then 
the event is placed in data array N(0,ip,E m i n ). The bin of 
zenith angle A9 = 1°, the bin of azimuth angle A<p = 0.5°, 
the bin of threshold energy E m [ n — 250 MeV. Integral muon 
intensity is calculated in the following way: 

T(R[nF , N(e,<p,E min ) 

1 ' £ SM ' £ MCS ' £add ' O Ji(0, <£, .fc m inj 

where N(9,tp,E m i n ) is the number of registered muons in 
a given angular and threshold energy bin. T is " live time" 
of registration. The parameter egM is efficiency of single 
SM triggering, and e ac id takes into account event rejection 
because of accompanying particles. Results of simulations 
and additional experimental data analysis give the follow- 
ing values: £sm = 0.936±0.004, e ac id varies from 0.83 to 



0.91 for different 9 and E m - ln (uncertainty of e a dd is less 
then 0.35%). The function SCl(0,(f,E m i n ) is the setup ac- 
ceptance calculated by means of MC method taking into 
account the structure of SM and selection requirements. 

Absolute muon intensity averaged in azimuth angle for 
zenith angles 61°<6><89° and for five threshold energies is 
represented in Table I and is shown in Fig. [3] (points) . Er- 
rors in the table include statistical and systematical uncer- 
tainties (uncertainty of threshold energy estimation, uncer- 
tainty of ssm, muon energy loss in the walls of surrounding 
buildings). 
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Fig. 3 

Dependence of absolute muon intensity on zenith angle at 
several threshold energies. symbols represent experimental 
results obtained in this work. curves are calculation 
results by formula (2). 

IV. Approximation formula 

For approximation of measured experimental data, the 
following simple formula is used: 

-^app (P l E m in ) — ~E^ 

■expf-^L-. E - -In (-MY ( 2 ) 
V 7+ 1 Ei+a-hx \hvcos6*JJ w 

The factor in front of the exponent reflects the form of 
muon spectrum in the upper atmosphere, and the expo- 
nential function takes into account muon decay. Here C\ 
is the normalization; E\ = E lran + a ■ (ho/ cos 9* — hi) is 
the threshold muon energy (GeV) at production level. In 
this formula a = 2.5T0~ 3 GeV-cm 2 /g is effective specific 
energy loss; (h / cos 9* — hi) is the path of muon in the 
atmosphere; ho — 1018 g/cm 2 is the total thickness of at- 
mosphere (altitude of setup under see level is taken into 
account); hi — 100 g/cm 2 is the effective depth of muon 
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TABLE I 

Dependence of integral muon intensity on zenith angle 



6° 




I{6)- 


10 5 , (cm 2 -s-sr) 








1.7 GeV 


3.3 GeV 


4.2 GeV 


5.4 GeV 


7.2 GeV 


61 


140.1 ± 2.3 










63 


122.5 ± 1.5 










65 


108.9 ± 1.2 










67 


93.6 ±0.9 










69 


79.9 ± 0.7 










71 


66.41 ± 0.6 










73 


54.8 ± 0.4 










75 


44.0 ±0.3 










77 


34.0 ±0.2 










78 


29.6 ±0.2 


25.5 ±0.4 








79 


25.4 ±0.2 


22.8 ±0.2 








80 


21.5 ± 0.2 


19.5±0.1 








81 


18.2 ± 0.1 


16.5±0.1 


15.7±0.4 






82 


14.9 ±0.1 


13.55 ±0.07 


12.9 ±0.2 






83 


11 98 ± 09 


10 98 ± 06 


10 78 ± 09 

1U. 1 U _L_ U.Vt/ 






84 


9.53 ±0.07 


8.79 ±0.04 


8.61 ±0.06 


8.1±0.1 




85 


7.27 ±0.06 


6.74 ±0.03 


6.56 ±0.04 


6.55 ±0.08 


5.84 ±0.09 


86 


5.44 ±0.12 


5.08 ±0.02 


4.97 ±0.03 


4.86 ±0.05 


4.59 ±0.04 


87 


3.95 ±0.09 


3.74 ±0.02 


3.65 ±0.02 


3.64 ±0.03 


3.51 ±0.02 


88 


2.78 ±0.06 


2.64 ±0.01 


2.61 ±0.02 


2.56 ±0.02 


2.50 ±0.01 


89 


2.01 ±0.05 


1.83 ±0.01 


1.74±0.01 


1.77±0.01 


1.69 ±0.01 



generation. E a = z -m/(c-T - cos 9*) is the effective crit- 
ical energy for muon; zq is the effective length at which 
the density of atmosphere is changed by a factor of e [13]; 
c is the velocity of light; t is muon life time and m is 
it's mass (GeV); cose?* = (cos a 6* + A") 1 '" is the approxi- 
mation of effect of atmosphere sphericity. As a result of 
fitting, the following values of free parameters were ob- 
tained: Ci = 0.103, 7 = 1.963, z = 6.81 km, A = 0.0577 
and a — 1.35. 

In Fig. 3, integral muon intensity calculated by formula 
(2) for five threshold energies (the curves) is compared 
with the present experimental data. Dependence of in- 
tegral muon intensity on zenith angle calculated for lower 
thresholds (1 and 0.3 GeV) and experimental data of ear- 
lier measurements [l]-[8] are presented in Fig. 4. Com- 
parison of calculated values with data [1,4,5,7,8] shows a 
reasonable agreement. In works [2] and [6], the intensity is 
somewhat higher than measured in [1] or calculated by (2). 
The integral intensity data at E m i n = 1 GeV obtained in [3] 
decrease with the increase of zenith angle more slowly than 
it follows from [1] and calculation by (2), but at angles less 
than 72° the agreement is quite well. 

V. Conclusions 

Experimental data of coordinate detector DECOR cover 
unexplored earlier region for integral muon intensity at 



threshold energies 1.7<-E m i n <7.2 GeV and zenith angles 
61°<6'<89 . It is important to mark that the measure- 
ments for all thresholds were performed simultaneously 
with a single setup, that minimizes systematic uncertain- 
ties. Extrapolation of the present data to lower thresholds 
is in a reasonable agreement with the results of other mea- 
surements. 
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Fig. 4 

Dependence of absolute muon intensity on zenith angle for 

LOWER THRESHOLDS (0.3 AND 1 GeV). CURVES ARE CALCULATION 
RESULTS BY FORMULA (2). SYMBOLS FOR EXPERIMENTS ARE THE 
SAME AS IN FlG. 1. 
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